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We analyse restrictions on the stochastic superspace parameter space arising from 1 fb _1 of LHC 
data, and bounds on sparticle masses, cold dark matter relic density and the branching ratio of the 
process B s — > /i + fi~ . A region of parameter space consistent with these limits is found where the 
stochasticity parameter, £, takes values in the range — 2200 GeV < £ < — 900 GeV, provided the 
cutoff scale is G>(10 18 ) GeV. 



Supersymmctry is a popular extension to the Standard Model, renowned for its solution to the hierarchy problem, 
while also providing a dynamical mechanism of electroweak symmetry breaking and a dark matter candidate in its 
minimal form (see 1] for a review). By postulating a symmetry relating fermionic and bosonic states, it predicts the 
existence of a superpartner of opposite spin statistics and degenerate mass for each Standard Model particle. These 
superpartners have not been observed, requiring that supersymmetry be softly broken for it to remain a viable model. 

A number of methods for breaking supersymmetry exist, which can be broadly categorised into two groups: those 
originating from a fundamental theory where supersymmetry is spontaneously broken in a hidden sector, then com- 
municated to the visible sector by a messenger, and the purely phenomenological approach, where no explanation is 
£NJ , provided as to the origins of supersymmetry breaking. Although fundamental theories provide a dynamical explana- 
tion for the mechanism of supersymmetry breaking, complications such as large sparticle induced FCNC amplitudes in 
gravity mediation, and CP violation issues in gauge mediation, leave no single compelling solution [2] . Alternatively, 
the phenomenological approach of writing down the most general Lagrangian with explicit soft-breaking terms fails 
to address the source of SUSY breaking, while introducing an overwhelmingly vast new parameter space. 

Stochastic superspace is a unique mechanism for softly breaking supersymmetry that bridges the two approaches 
discussed above Q. By considering the Grassmannian coordinates to be stochastic variables, a very constrained 
set of soft-breaking terms emerges. Though the underlying cause of stochasticity is not postulated, it introduces 
predictability not achievable through the purely phenomenological approach. 
Qh| In this paper we investigate the phenomenology of stochastic superspace at the LHC. For each suitable point in the 
stochastic superspace parameter space, we use a fast simulation of the ATLAS detector to investigate whether points 
i ' are excluded by the ATLAS zero lepton searches, which in turn are amongst the most constraining current limits on 
direct sparticle production. 

The rest of the paper is structured as follows: Sec. Uprovides an overview of stochastic superspace, Sec. HIl discusses 
constraints applied in a scan of its parameter space, and Sec. IIIII describes the technique used to simulate stochastic 
superspace model points at the ATLAS detector, and explains how we approximate the ATLAS exclusion limits on 
\ supersymmetric particle production. We terminate with concluding remarks in Sec. IIVI 



(N 



43 ■ 

Or 



o 

(N 



X 



I. OVERVIEW OF STOCHASTIC SUPERSPACE 

The basis of stochastic superspace models is the assumption that the Grassmannian coordinates, and 9, are 
stochastic variables. Writing down the most general probability distribution consistent with Lorentz invariance, we 
find that only one additional parameter is required such that [3j, 

V (6, 9) |C| 2 = 1 + T (99) + £ (06) + (99) (90) , (1.1) 

where £, a complex number of mass dimension, is the measure of stochasticity. The Lagrangian in ordinary spacetime 
is then found by averaging the supersymmetric Lagrangian over the probability distribution in Eq. (11.11) . Applying 
this procedure to the superpotential of the minimal supersymmetric standard model, 

Wsm = nH u H d + y up QU c H u + y dov > n QD c H d + y lcpt LE c H d , (1.2) 
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leads to the soft-breaking terms 



L so ft scalar — ~ £* I^H u Hd — 2£* 



y up QU c H u + y down QD c H d + y lcpt LE c H d 



+ h.c, (1.3) 



where tildes represent the scalar component of the chiral superfield. Similarly, averaging over the gauge-kinetic F 
densities results in soft-breaking masses for the gauginos, x , 



igaugc 



5-^TrA w A W +h.c, (1.4) 



where W^ a denotes the standard model field-strength superfields. See Q for further details. Referring to Eqns. 
(|1.3p and (|1.4[) . it is clear that the tree- level soft-breaking terms for minimal stochastic superspace are [3j: 

(i) the bilinear Higgs soft term, = £*, 

(ii) the universal trilinear scalar soft terms proportional to the Yukawa couplings, Aq = 2£*, 

(iii) the universal gaugino masses, mi/2 = \ 

(iv) the universal soft scalar masses, TOq = 0- 

As these soft-breaking terms are renormalised at the quantum level, the relations listed above are only defined as such 
at some energy cutoff scale, A, which becomes an additional parameter of the theory. 

The emergence of a pattern of soft-breaking terms dependent only on the parameter £ categorises stochastic su- 
perspace models as a special case of the broader constrained supersymmetric standard model (CMSSM). Although 
models with absent soft scalar masses are typically considered excluded 1 , this is because the soft-breaking parameters 
are customarily defined at A = Mgtjt ~ 10 16 GeV. We have found that for larger cutoff scales, A > 10 18 GeV, the 
neutralino becomes the LSP, and consequently, one uncovers a new phenomenologically viable region of parameter 
space We look at this in further detail in Sec. HI] by computing the neutralino relic density. 



II. SCAN OF PARAMETER SPACE 



We have previously demonstrated the existence of phenomenologically viable regions of the (£, A) parameter space 
falling within the stau coannihilation regime 0, Q • In this study we expand the scope of the analysis primarily to 
examine a broader range of £ values and include a calculation of the cold dark matter relic density. 

Most publicly available sparticle spectrum software packages take the mSUGRA parameters mo, m 1 / 2 , Ao, tan/3 
and sgn/x as inputs at the cutoff scale. However, stochastic superspace makes a prediction for the bilinear Higgs soft 
term j3 M , rather than explicitly predicting a value for tan/3. To facilitate a more thorough analysis of the model, 
we previously modified the program SDFTSUSY 3.1.7 [f| to swap tan /3 for as an input. As branching ratios are 
required for this analysis, we use the package Isajet 7.81 [7[ to calculate the sparticle and decay spectra, using the 
value of tan/3 found with SOFTSUSY. Various packages included in IsaTools calculate the relic density and branching 
ratios of rare decays. 

We scan the parameter regions 

-2600 GeV < £ < -200 GeV, (2.1) 
M GU T < A < Mm, (2.2) 



1 This is because under the assumption of R-parity conservation the lightest stau is a stable lightest supersymmetric particle (LSP) with 
the mass below the CDF model-independent limit on the mass of a charged massive stable particle (CHAMP), moHAMP ^ 250 GeV at 
95% C.L. [3 
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FIG. 1: The phenomenologically viable region of parameter space, where the black area shows points that satisfy the WMAP 
bounds on relic density, and the gray area shows points where the relic density is lower than observed, but satisfies all other 
constraints. 



over 150 points in £ and 50 points in A for a total of 7500 points, against the following conditions: 



m^a < m-fj, (2.3) 

m h >109GeV, (2.4) 

m ffl > 95.7 GeV, (2.5) 

m~ bR > 89 GeV, (2.6) 

m? R > 81.9 GeV, (2.7) 

BR(B S -> fi+fi') < 1.2 x 10" 8 , (2.8) 

0.0988 < VLh 2 < 0.1252. (2.9) 



In stochastic superspace, the character of the LSP is dependent on one's parameter choice, and will either be the 
lightest neutralino or the lightest stau. Since R-parity is assumed to be conserved, we enforce the condition in Eq. 
(|2.3p to ensure all instances of a stau LSP are excluded. The lightest CP— even Higgs must satisfy condition (12.41) . 
where we take into consideration the theoretical error of 3 — 5 GeV 8] in the published limit of ?>?/, > 114.4 GeV at 95% 
CL Q. Conditions (|2.5[) - (I2.7[) arise from supersymmetric particle searches [9|. Condition (|2.8[) . the upper limit on the 
branching fraction of the process B s — > is as measured by the CMS and LHCb collaborations [l(| El- Finally, 

condition (|2.9[) shows the acceptable range of relic density based on the WMAP 7-year mean, flh 2 = 0.1120 ± 0.0056 
[lH, and the theoretical error of ±0.012sijsy 0, shown with errors added in quadrature. 

Applying all constraints listed in Eqns. (I2.3I) - (|2.9[) . we find a phenomenologically viable region of parameter space, 
shown in black in Fig. \T\ Relaxing condition (|2.9[) . we find another region of interest where Vth 2 < 0.0988 (see 
gray region in Fig. [1]). Although the relic density is lower than observed, this region is not necessarily excluded 
as the minimal model of stochastic superspace may be extended to include an additional source of dark matter. 
Furthermore, theoretical uncertainties in the determination of the sparticle mass spectrum can have a significant 
effect on the calculated value of the relic density, where a 1% difference in mass could affect the calculated relic 
density by up to 10% [H], [HJ- It is possible, therefore, that the region of parameter space satisfying the relic 
density experimental bounds may be larger than depicted in Fig[TJ The region where A < 1.15 x 10 18 GeV and/or 
£ > —600 GeV is excluded due to sparticle masses violating conditions (I2.5l) - (|2.7p or going tachyonic. Excluded points 
in the region where £ < —1200 GeV and A > 1.15 x 10 18 GeV have relic density larger than the upper bound in 
condition (|2~9"|) . In the region with -1200 GeV < £, < -640 GeV and A > 1.15 x 10 18 GeV, a number of points arc 
excluded for having a stau LSP. For clarity, Figure [2] displays these regions graphically. All points in the parameter 
space satisfying constraints (|2.3p - (|2.8[) above, with relic density within or below the bound given in condition (|2.9p . 
form the phenomenologically viable region of stochastic superspace. 
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FIG. 2: The causes of data point exclusion in the parameter space scan is shown in this plot. The viable region shown in 
gray is surrounded by three types of exclusion zones: too small or tachyonic sparticle masses are responsible for the exclusion 
of points in the region labelled "sparticle masses," relic density higher than the upper WMAP bound excludes points in the 
region labelled "relic density," and a stau LSP is responsible for excluding points in the region labelled "LSP." 



III. IMPACT OF LHC EXCLUSION LIMITS 
A. LHC SUSY searches 

The ATLAS and CMS experiments [l6l [ljj have recently updated their searches for supersymmetric particles 
using the 2011 dataset [LSl - l25l ]. Data collected from proton collisions at the Large Hadron Collider at yfs = 7 TeV are 
analysed in a variety of final states. No excess over the standard model expectation has been observed, allowing one to 
set exclusion limits in a variety of candidate model spaces such as the CMSSM, the model that most closely resembles 
that considered here. As the LHC is a proton-proton collider, one expects to dominantly produce coloured objects 
such as squarks and gluinos, whose inclusive leptonic branching ratios are relatively small, and hence the strongest 
CMSSM exclusions to date result from the ATLAS searches for events with no leptons and the CMS searches for 
sparticle production in hadronic final states. The ATLAS and CMS limits have a similar reach in the squark and 
gluino masses, and here we consider only the ATLAS zero lepton limits for simplicity. Recent interpretations of LHC 
limit results can be found in [l^, I26l - l30j . 

ATLAS defined a series of signal regions in which to look for sparticle production, each tuned to enhance sensitivity 
in a particular region of the TOo-ttxx/2 plane. Events with an electron or muon with px > 20 GeV were rejected. 
Table U summarises the remaining selection cuts for each region, whilst Table |H] gives the observed and expected 
numbers of events. These numbers were used by the ATLAS collaboration to derive limits on a x A x e, where a is the 
cross-section for new physics processes for which the ATLAS detector has an acceptance A and a detector efficiency 
of e. These results are also quoted in Table [TT1 

The ATLAS collaboration has used the absence of evidence of sparticle production in 1 fb _1 of data to place an 
exclusion limit at the 95% confidence level in the too- to i/2 plane of the CMSSM for fixed Aq and tan/3, and for fi > 0. 
Although the stochastic SUSY model considered here may be considered a subset of the CMSSM, it is non-trivial 
to recast this limit into a constraint on the parameters £ and A. Given a signal expectation for a particular model, 
however, one can easily evaluate the likelihood of that model using the published ATLAS background expectation 
and observed event yield in each search channel. By simulating points in the £-A plane, we can therefore investigate 
the LHC exclusion reach in stochastic SUSY space, provided that we can demonstrate that our simulation provides 
an adequate description of the ATLAS detector. 

B. Simulation details and validation 

In this paper, we use Isajet 7.81 Q to produce SUSY mass and decay spectra then use Herwig++ 2.5.1 [3l[ 
to generate 10,000 Monte Carlo events. Delphes 1.9 [32[ is subsequently used to provide a fast simulation of the 
ATLAS detector. The total SUSY production cross-section is calculated at next-to-leading order using PR0SPIN0 2 . 1 
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Region 


Rl 


R2 


R3 


R4 


RHM 


Number of jets 


> 2 


> 3 


> 4 


> 4 


> 4 


Ef ss (GeV) 


> 130 


> 130 


> 130 


> 130 


> 130 


Leading jet pr (GeV) 


> 130 


> 130 


> 130 


> 130 


> 130 


Second jet pr (GeV) 


> 40 


> 40 


> 40 


> 40 


> 80 


Third jet p T (GeV) 




> 40 


> 40 


> 40 


> 80 


Fourth jet pr (GeV) 






> 40 


> 40 


> 80 


A<Mjet, pT ss Uu 


> 0.3 


> 0.25 


> 0.25 


> 0.25 


> 0.2 


m e fi(GeV) 


> 1000 


> 1000 


> 500 


> 1000 


> 1100 



TABLE I: Selection cuts for the five ATLAS zero lepton signal regions. A0(jet, p™ lss ) m in is the smallest of the azimuthal 
separations between the missing momentum p™ lss and the momenta of jets with pr > 40 GeV (up to a maximum of three in 
descending pr order). The effective mass m c g is the scalar sum of _E™ 1SS and the magnitudes of the transverse momenta of the 
two, three and four highest pr jets depending on the signal region. In the region RHM, all jets with pr >40 GeV are used to 
define m e g. 



Region 


Rl 


R2 


R3 


R4 


RHM 


Observed 
Background 
a x Ax e (fb) 


58 

62.4 ±4.4 ±9.3 
22 


59 

54.9 ±3.9 ±7.1 
25 


1118 
1015 ± 41 ± 144 
429 


40 

33.9 ± 2.9 ± 6.2 
27 


18 

13.1 ± 1.9 ± 2.5 
17 



TABLE II: Expected background yields and observed signal yields from the ATLAS zero lepton search using 1 fb 1 of data [T^ ]. 
The final row shows the ATLAS limits on the product of the cross-section, acceptance and efficiency for new physics processes. 

[33j . where we include all processes except direct production of neutralinos, charginos and sleptons since the latter 
are subdominant. The ATLAS setup differs from this only in the final step of detector simulation, where a full, GEANT 
4 based simulation 34] is used to provide a very detailed description of particle interactions in the ATLAS detector 
at vast computational expense. 

It is clear that the Delphes simulation will not reproduce every result of the advanced simulation. Nevertheless, 
one can assess the adequacy of our approximate results by trying to reproduce the ATLAS CMSSM exclusion limits. 
We have generated a grid of points in the mo-mi/2 plane using the same fixed values of tan/3=10 and Aq=0 as the 
published ATLAS result. We must now choose a procedure to approximate the ATLAS limit setting procedure. 
ATLAS use both CL S and profile likelihood methods to obtain a 95% confidence limit, using a full knowledge of the 
systematic errors on signal and background. Although the systematic error on the background is provided in the 
ATLAS paper, we do not have full knowledge of the systematics on the signal expectation, which may in general 
vary from point to point in the mo-m 1( / 2 plane. Rather than implement these statistical techniques, we take a similar 
approach to that used in [3(|, and use the published a x A x e limits to determine whether a given model point is 
excluded in a search channel. We use our simulation to obtain the a x Ax e value for a given model point, and consider 
the model to be excluded if the value lies above the limit given in Table [TT1 This allows us to draw an exclusion contour 
in each search channel, and we estimate the combined limit by taking the union of the individual exclusion contours 
for each channel (i.e. the most stringent search channel for a given model is used to determine whether it is excluded) . 

The procedure defined above neglects systematic errors on the signal and background yields and, as noted inj30j, 
this leads to a discrepancy between the Delphes results and the ATLAS limits in each channel. We follow [30| in 
using a channel dependent scaling to tune the Delphes output so that the limits in each channel match as closely 
as possible "by eye". We obtain factors of 0.82, 0.85, 1.25, 1.0 and 0.70 for the Rl, R2, R3, R4 and RHM regions 
respectively. Comparisons between the resulting Delphes exclusion limit and the ATLAS limit are shown in Figure|3l 
where we observe generally good agreement in all channels. The largest discrepancy is observed in the RHM channel, 
where we find that one cannot get the tail of the limit at large uiq to agree with the ATLAS limit whilst simultaneously 
guaranteeing good agreement at low mo. This is likely to be due to the fact that we have effectively assumed a flat 
systematic error over the mo — mi/2 plane, whereas the ATLAS results use a full calculation of the systematic errors 
for each signal point. It is important to notice however that the combined limit will be dominated by regions Rl and 
R2 at low mo, and thus by choosing to tune the RHM results in order to reproduce the large mo tail, one can ensure 
reasonable agreement of the combined limit over the entire range. Where disagreement remains, the Delphes limit 
is less stringent than the ATLAS limit, and hence using it gives us a conservative estimate of the ATLAS exclusion 
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FIG. 3: Comparison between Delphes and ATLAS 95% exclusion limits in the mo — TO1/2 plane, for the signal regions Rl, R2, 
R3, R4 and RHM defined in Table [I] In the combined limit plot, the ATLAS limit is obtained using the ATLAS statistical 
combination, whilst the Delphes limit is obtained by taking the union of the Delphes limits for each signal region. 



reach. 



C. Limits on stochastic SUSY space 



Having obtained reasonable agreement with the ATLAS GMSSM limits in the previous section, we now turn our 
attention to calculating the LHC exclusion zone in the plane of stochastic SUSY parameters £ and A. To avoid 
unnecessary (and computationally expensive) simulation, we only simulate those points that passed the constraints 
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FIG. 4: ax Axe values for stochastic SUSY points, as obtained using the simulation code detailed in the text, for ATLAS signal 
regions Rl, R2, R3 and R4. The acceptable region of stochastic supersymmetry parameter space can be found by referring to 
the ATLAS limits shown in Table [LTl Points on the graphs above with values exceeding those limits are excluded. For example, 
the limit for signal region Rl is 22 fb. All stochastic superspace points with a x A x e larger than 22 fb are ruled out, in this 
case £ ^ —850 GeV, represented by the region to the right of the light blue colour. Correspondingly, the limit for R2 is 25 fb, 
excluding £ > -800 GeV in the region to the right of the light blue colour. The limit for R3 is 429 fb, which does not exclude 
any stochastic superspace points. R4 has a limit of 27 fb, which excludes £ > —750 GeV in the region to the right of the green 
coloured squares. 



detailed in Sec. ILTl 

Figures S] and [S] show the ax Ax e values (including the systematic factors) for each of the ATLAS search channels. 
a x Ax e is most strongly dependent on £ which is to be expected; points at low |£| have light coloured sparticles and 
will thus have large production cross-sections at the LHC. Raising the mass scale essentially reduces the production 
cross-section leading to a corresponding decrease in the a x A x e value as one travels left along the £ axis. There 
is little to no dependence of a x A x e on the cutoff scale for a given value of £. We can conservatively bound the 
region £ < —900 GeV as consistent with 1 fb -1 of LHC data based on ATLAS signal region Rl, which is the most 
constraining in the case of stochastic superspace. 

IV. CONCLUSION 

In this paper, we perform a systematic scan of the parameter space of the stochastic superspace model of soft 
supersymmetry breaking. Its phenomenology is checked against constraints from direct sparticle searches, cold dark 
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FIG. 5: a x A x e values for stochastic SUSY points for ATLAS signal region RHM. The ATLAS limit for this signal region is 
17 fb, which excludes stochastic superspace points with £ > —700 GeV in the region to the right of the light green colour. 



matter relic density and the branching ratio of the process B s — > fi + /i~ . Points that pass these conditions are analysed 
further, obtaining a value for a x Ax e, which can be directly compared to ATLAS search limits from the zero lepton 
channels. We find that the stochasticity parameter is restricted to the region £ < —900 GeV, while the viable region 
of the cutoff scale, A, shrinks with increased |£|, bounded above by relic density concerns, and below from sparticlc 
mass constraints. 

Since a prime motivation for the the low-energy supersymmetry is a solution to the hierarchy problem, the increased 
bounds on sparticle masses seriously undermine its validity. The required tuning of parameters is estimated at the 
level of ~ 10% for £ ~ —1000 GeV @. Thus, stochastic supersymmetry still allows a region of parameters that 
provide a satisfactory solution to the hierarchy problem. However, as has been shown in 0, the mass of the lightest 
Standard Model-like Higgs boson is predicted in a very narrow range ps 112 — 116 GeV in the minimal model. 
This range is outside the 124 — 126 GeV region for the Higgs mass hinted in the very recent ATLAS [35[ and CMS [36[ 
data. If the latest evidence for the Higgs mass will be confirmed, the minimal model with stochastic supersymmetry 
will certainly be excluded, motivating research into a non-minimal implementation of stochastic superspace consistent 
with these results. 
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